Intermediate spin, Schrodinger cat states and nano- magnets. 
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Quantum tunneling of nano-magnets finds a natural description in terms of intermediate spin. 
Periodic magnetic effects correspond to a change of flux by the flux quantum $o. Schrodinger cat 
states with different superpositions of the applied magnetic fleld occur. The molecular magnet Feg 
is discussed and oscillations are predicted for Mni2. 
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Recent experiments by WernsdofFer and Sessoli fit] on 
the molecular magnet known as Feg show fasinating be- 
haviour as a function of an applied field. This is an or- 
thrombic easy axis material with one perpedicular axis 
harder than the other. As a function of either an easy or 
hard axis field the tunnel splitting manifests self-similar 
behavior at periodic values of these fields. As a func- 
tion of the third component the splitting increases es- 
sentially exponentially. While these results can be well 
reproduced by the exact diaginalisation of a suitable large 
spin model, as yet there is not an analytic theory which 
permits a clear understanding of the observed phenom- 
ina. An adequate such theory should permit the predic- 
tion of new related effects for materials, e.g., of different 
symmetry. 

Invariably, in the tunnelling regime, the least action 
paths which make any significant contribution to the tun- 
neling amplitude can be mapped to a circle. The problem 
is then equivalent to the motion of a charged particle in 
a magetic field. Fig. (Q), i.e., to the problem of interme- 
diate spin. In this Letter will be described an analytic 
theory based upon an extension of the intermediate spin 
approach developed earlier by the author H. In terms 
of this the behavior of Feg is easily understood and, e.g., 
a new periodic effect can be predicted for the large spin 
model often used to describe the much studied |^ Mni2. 

The relevant spin model for Feg can written as: 



H = -iD-E)Sj + 2ES^'^ 

—ES — hSz — hiSx — htSy, 



^/^o = h/4:hy/2E{D + E) evident in experiment |]l|. 

There are some interesting surprises, e.g., in order to 
describe tunneling in Feg it is necessary to construct 
Schrodinger cat states which involve different values of 
the applied magnetic field. Even more bizarre is that the 
other transverse field ht is equivalent to an imaginary 
fiux $ leading to a cosh dependence of the splitting. In 
general for an arbitrarly directed transverse field the ef- 
fective fiux is complex! 

There are classical ground states which with h — ht — 
lie along the positive and negative x-axes. Tunneling 
occurs between these ground states. Adding a the longi- 
tudinal field hf^, e.g., the ground state for one direction 
can be brought to the same energy as the first excited 
state in the other direction. In this case one wave func- 
tion has even parity and the other odd and this adds a 
phase shift of tt to the wave function. This also con- 
verts whole into half-integer spin, again an effect seen in 
experiment [nj. 

These ideas permit the prediction of new effects even 
for the extensively studied molecular magnet Mni2. This 
is a tetragonal easy y-axis magnet but with equivalent x 
and y hard directions in the transverse plane. There are 
four equivalent least action tunnel paths in the absence 
of a field. When a modest field is applied along any of 
the four hard directions there remain only two equivalent 
least action paths which are mapped to the intermediate 
spin circle. Actually this problem can be solved even for 



(1) 



where the positive anisotropy constants have E < D 
iQ. The magnetic field h, parallel to the hardest axis, 
is equivalent to that in Fig. (II|) and most physical quan- 
tities, including the tunnel splitting, must be periodic in 
the flux $ produced by this field, this with a period $0i 
the flux quantum. If the physical spin value S is whole 
integer, a flux of {n + ^)^o corresponds to a statistical 
parameter ag ~ 1, modulo unity, and transmutes the 
model to half-integer spin while intermediate values n$o 
reflect a a^ = 0, modulo unity, and are whole-integer 
points. Kramer's theorem implies a zero tunnel splitting 
for half-integer points. These facts are reflected in the 
quasi-periodic tunnel splitting A = Aq cos 27r$/$o with 




FIG. 1. A particle of charge q is confined to a unit circle. 
A tunneling problem might be defined by adding a potential 
V such that V{<l)) ~ V{4> + tt). If this potential has minima 
at = and n and the barrier between is higher than the 
zero point energy then there will be a ground doublet with 
a splitting A. A perpendicular field exerts no force on the 
particle but does change the slitting as explained in the text. 
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FIG. 2. The analytic result is Ao| cosli(7r/i/p) cos(7r/i/p) 
where p ~ 2S^'^(12aD^)^'''. This corresponds to a complex 
magnetic field in which the real and imaginary parts are equal. 
This exact result is indicated. The difference becomes smaller 
with increasing S. 

small fields by mapping the four paths to two circles and 
then adding the wavefunctions in much the same way as 
in the elementary theory of a two slit experiment. These 
circles have the applied field at 45° and the flux is com- 
plex with equal real and imaginary parts. As with Fcg 
there are periodic half-integer points at which the split- 
ting is zero but now the maximum splitting increases 
by almost an order of magnitude per period. The re- 
sult is given in the caption to Fig. (0) which compares 
it with that of exact diagonalisation. Further details will 
be given elsewhere [0. It can also be shown that similar 
oscillations occur when the field is directed along a hard 
three or four fold axis of a cubic material B . 

For the basic intermediate spin problem H of a charged 
particle in two dimensions, Fig. (ffl), the relevant group is 
S0(2) which is equivalent to U(l), i.e., is abelian with a 
single operator 5*2, the generator of rotations e*"^^"^, where 
the z-axis is perpendicular to the circle and where (f) is 
the angle to the a;-axis. For a unit circle Sz = p,f, the mo- 
mentum conjugate to 0, i.e., [0,P0] = ih. The eigenstates 



of Sz and p^ are ■(/'m(</>) 



^im<f) 



but for particles on 



(27r)i/2< 

a circle there is no reason to insist on m being an integer 
reflecting the fact that the spectrum of Sz for 5*0(2) is 
continuous [g. In general for this singular gauge 'ipm{4>) is 
multi-valued. A convenient single valued function can be 
defined by taking some convention for ipm {4') in the inter- 
val {0, 27r} and then analytically continuing this function 
to the interval {— oo,+oo}. For any physical potential 
V{4)) = V{4> + 2?™); n = 0, ±1, ±2 . . ., and implies that 
only the basis states ';/'2^+n(0); "■ = 0,±1,±2..., mix. 
It follows for a given statistical parameter as the energy 
spectrum is discrete. This can be stated differently: Since 
V{4)) is periodic the solutions obey Floquet's (Bloch's) 
theorem, i.e., V'fc(</') = e*'^*Mfc((/)) where Uk{(j)) is periodic. 
The reciprocal lattice vector K — 2n/{a = 2tt) = 1 and 
the energy e(fc) = €{k + nK); n — 0,±1,±2... i.e., 
is periodic. An examination of the Fourier components 
immediately implies k = as/2. Given that V{(j)) has a 
point with reflection symmetry the ip±ki4') ^-re degener- 
ate and the most general solution for a given energy is 




FIG. 3. Mapping of the unit sphere to a two dimensional 
rectangle 

Aipk{4')+Bip_i.{<j)) for which the spectrum of m — n±-^, 
i.e., there are twice the number of observable Sz values 
except importantly for the whole and half-integer points 
when Qfg = or 1. This more general method of extrapo- 
lating between whole and half-integer spin which will be 
found applicable here. 

In the non-singular gauge the basis set V'n(^) ^iels 
n = ±1, ±2 . . . and the charged particle encircles a mag- 
netic flux $ = as$o/2 where $o is the appropriate flux 
quantum. For the same V{(l)) there exists an eigenstate 
Uk{(j)) with the same energy as in the singular gauge, 
i.e., the energy spectrum are identical and two prob- 
lems can mapped to each other and both reflect inter- 
mediate spin. In the non-singular gauge the combina- 
tion Auk{<j)) + Bu^kiff)) comprises a linear combination 
of states with the flux <I> in both its positive and negative 
senses. Evidently the expectation value of the flux ($) 
is in general smaller than as^o/2. This possibility of a 
linear combination of different flux states does not seem 
to have been envisaged in connection with intermediate 
statistics although it would arise if, e.g., quantum coher- 
ence in SQUIDS could be observed. It is realized in the 
systems, e.g., Fes, discussed here. 

The usual functional integral formulation n involves 
coherent state \9,(j)) obtained by rotating \Sz = S) 
through usual angles 9,4>. However the |0, 0) are very 
much over-complete. In fact the states |0) = \0 = f , </>) = 
{2S + l)^/'^Y.l^^_se-'""t'\m), where \m) = \Sz = m), 
remain over-complete. If \^) = X]m=-s '^('™)l'™) then 



^((/.) = ((/-iv) = e; 



^ira<l} 



a{m), i.e., is a Fourier 



transform of the TTi-space wave function a{m). The 
function ip{4>) is defined on the interval {— cxd,-|-cx3} 
and the inverse transform ip{m) = 2tt I~oo '^4>'4'{4>) — 
X^m' <5(m — m')a(m') so that |V') = /_oq rfm-0(m)|m), 
i.e., ij}{m) is a definition of the m-space wave function 
appropriate when m is considered as a continuous vari- 
able. Clearly Szip{(f>) = J_ dme'^™''^hmip{m) — p^ip^cj)) 
where p^ = —ih-jx is the same as that defined ear- 
lier. Also evident is that multiplication by e*"^ amounts 
to a displacement of m by unity and so, e.g., S""^ = 
e^'''[h^S{S+l)~P4,{p4, + 1)^1^. These definitions permit 



an arbitrary spin Hamiltonianto be expressed in terms of 
the cononical p^ and (p. 

In terms of the unit sphere (j) and p^ might be consid- 
ered as a particular projection in which the axis of quan- 
tization plays the role of the North-South axis as shown 
in Fig. (ph. In these terms, the angle is evidently the 
longitude while p^ is the projection on the z-axis and is 
therefore the sine of the latitude. 

It is necessary to carefully formulate the problem in 
terms of the basis set \m), i.e., in ?Ti-space. If IV') = 
^^^_g a(?7i)|?7i), then in this basis Schrodinger's equa- 
tion for Eqn. (||) can be written as {hi — ht = 0): 

(e - {2Eh'^m'^ - hhm)) a{m) 






-iAC:^(m + 2) 



+ [M, 



m+l2 



MI: 



i{m - 2) 
'^ ]a{m) 



(2) 



where M* = [S{S + 1) - st]^^^. The structure of this 
equation is worth noting. It has the form of two tight 
binding models, one for even m values and the other for 
odd m. 

The wave function for the eigenstate \Sx — S) is 

nl/2 



a'^im) 



W- 



(2g)! 



and a'^{m) has to be a 

E — 0, this corresponding 

Using Stirling's approxi- 

however the 



solution to Eqn. (0) with h = 
to the absence of tunneling 
mation this reduces to a^{m) oc e~™ /^'^ 
deviations from this approximation are very important. 
The absence of tunneling when h = E = is put in ev- 
idence by a transformation: aim) = a^{ra)f{ra). With 
the definition /(0) = ^^ 
Schrodinger's equation: 



e*'"'^/(m) this yields the exact 



(£-2i?p^2y(0) = y(0,p^)/(0) 



(3) 



with 



V(<i,,p^)^\{D-E)y^ 
2 sin^ (j){h^S{S -l)+p^^)- ih{2S - 1) sin 2(/)j5 J , (4) 



where f = {e-E°); E° ^ ~h^S^{D + E)). This is of the 
intermediate spin form but with a momentum dependent 
potential, V{(j),Pci,). 

A new wave function is defined by /(0) = (1 ^ 
{asP(j,/2hS))g{4)). This admixes, e.g., a little of the first 
excited into the ground state. If g{(l)) is the solution of 

(e - 2E {p^ + h^f WE^)g{cj,) 

^V{c^,p^ + h^)g{c^), (5) 

then for modestly large S, f{4>) = (1 — {asp^/2hS))g{(l)) 
is the solution to 



(e - 2Ep/ 



2h{D 

= V{(f>,p^)f{^). 



E)YP^]fi4') 



(6) 



The shift in p^ determines the wave vector: k — as/2. 
In turn the statistical parameter must be chosen to 
give ground state energy determined without tunnel- 



ing. Required is that h?E^ — 



8(D+E) ' 



so fc = 



^ = h/2h^2E{D + E). This fixes the Zeeman term 
in Eqn. (|); 2h{D + E)^p^ = y/{D + E)/2Ehp4,. 
Since, in general, this is larger than hp^ the solution 
is a Schrodinger cat state, A{1 — {asP^/2hS))uk{(j)) + 
B{1 + {asP,p/2hS))u-k{4>) involving a linear combina- 
tion of both bare fields ± ^/{D + E)/2Eh such that the 
expectation value is equal the physical applied field. 
This requires A = 1 + ^J2E/{D + E) and B = 1 - 
\J2E lip + E). The total energy for level n and wave 
vector fc is. 



S„(fc) = S° 



^^^^-^+^n-H-);..o + e„(fc). (7) 



The first two terms are purely classical in origin while 
the third corresponds to the harmonic motion about the 
bottom of the well. It is the last term which reflects the 
tunneling between wells. 

The problem therefore reduces to the determination of 
e„(fc) and when the tunneling approximation is well jus- 
tified this reflects the tight binding approximation. With 
only a transverse field /i, the problem has a higher sym- 
metry and potential y(0,p0) has a ^ period of a = tt 
rather than 2-k so the reciprocal space unit vector K ^ 2. 
For a tight binding model the tunneling energy, e„(fc) = 
^ cos Trfc where Ao/2 is the undetermined matrix ele- 
ment for tunneling for level n. (The method for deter- 
mining this is given in y.) There are two solutions which 
are acceptable for h = 0, namely fc = ^ = corre- 
sponding to even m values and fc = ^ = 1 and odd to. 
These two k values correspond to e(fc) — ±-^ and so 
the zero field tunnel splitting is reflected in the param- 
eter Aq. For finite transverse fields h the correspond- 
ing fc = % 



'f = {h/2hy/2E{D + E)) and fc = 1 -h ^ 
1 + (h/2hy^2E{D + E)) and so the tunnel splitting is. 



A = Aq cos 



nh 



2hy/2E{D + E) 



(8) 



i.e., of the form mentioned in the introductory remarks. 
An equivalent result was first obtained some time ago by 
Garg and attributed to "topological quenching" . 

Nowhere in the development has use been made of the 
assumption that 5 is a whole- integer. The only difference 
when S* is a half-integer is that the values of m are also 
half-integer and thus the solutions correspond to fc = ^ -|- 



^ = i + {h/2hy/2E{D + E)) and fc = -i + ^ = -i + 
{h/2hyj2E{D + E)). This has the effect of replacing the 
cosine is replaced by a sine. 



Consider first the effect of adding hi alone. This adds 
a potential, i.e., V{4>iP<p) =^ ^(^iP^) + fiShicoscj) and 
removes the symmetry between cf) — Q and tt. The 
period of the potential is now the full a = 27r, with 
K = 1, hut with two wells per unit cell. The har- 
monic levels near <j) = tt have quantum numbers des- 
ignated by n and are higher than those near ^ = 
with labels n' . Ignoring tunneling these have energies. 



hnhf 



En - Eo - s(D+E) 
^0 - s{D+E) - ^'^'^t 

is only one solution with k 



_ + (n + ■^) fixjQ, and, £"„/ = 

+ (n' + i) fixjQ. For a given h there 

^ = {h/2h^2E{D + E)). 
There are now two possibilities (i) £"„ ^ En' and there 
is a very narrow band formed about both these energies. 
The energy differences have a very small dependence on 
k and hence /i, see below, (ii) When En ~ En' there 
is resonant tunneling. There are again two cases (iia) 
when both n and n' are odd or even and (lib) when one 
of n and n' is odd and the other even. As illustrated 
in Fig. (W) there are two tunneling matrix elements one 



for the barrier between 



and TT and A^ 



for 



that between cj) = —tt and (or tt and 2tt). With first 
h — {ht =)0, Schrodinger's equation admits real solu- 
tions, 'ip{m) or V'(0)j s-iid so these matrix elements can 
be made real. For case (iia) when the symmetry of the 
wave functions is the same the tunneling matrix elements 
A„^„' = A„/_„ = Aq. For finite ft,, for the Floquet (Bloch) 
wave function, the phase advances by 2TTk from one cell 
to the next while to agree with the smaller unit cell ap- 
propriate when hg — the phase should advance by tt/c 
from one well to the next within the cell. This leads to 



A. 



e'^'^'An and A„ 



-iTrk 



Ar 



At this point it is possible to introduce the second 
transverse field ht by slight of hand. It is observed 
that if E changes sign then the two transverse axes 
change role and thus h becomes the equivalent of ht- 
If H is considered to be a function of a complex pa- 
rameter E the solution might be analytically contin- 
ued from positive to negative values. In particular 
k 
k 



{h/2hy^2E{D + E)) becomes pure imaginary, i.e., 
i{ht/2hy/2E{D 



ikf — i[ 

= e^^ ('=+*'=' )Ao with A^ 



E)) and it follows that 
„ = e-*^('=+''=*)Ao. This 



information can be represented by a two level model: 



n=iEn-En')S, 



Ao r 



^i-n{k+ikt) c+ i -i-n{k+ikt) o- 



where the center of gravity energy ^ (£'„ -I- En' ) has been 
dropped. Degenerate perturbation theory then yields. 



ek ^±\J{En - En'Y + AAQ^\cosTT{k + ikt)Y or at res- 
onance €k = ±2 I cos7r(fc -I- ikt)\, with k 

and kt — 



2 

ht 



2hy/2E{D+E) 

The hi = ht = result is recov- 



2hy/2E{D-E) ■ 

ered when n = n'. When for (lib) one wave function 




FIG. 4. When, e.g., n — and n' — 2, both wave functions 
are even there is no sign change. If n = and n' = 1 the wave 
must change sign as its (j) argument advances by 27r. 

changes sign while the other does not it is necessarily the 
case, for h = ht = 0, that A „/ — — A„/ „ = Aq whence 
everywhere in ej. the cosine becomes a sine and the result 
is as if the original S was half-integer! This explains the 
alternation observed in experiment m. 

Certainly the most striking claim is that the wave iunc- 

tioiiA{l-{asP4,/2hS))uki4>) + B{l + {asP4,/2hS))u^ki4>) 
for an intermediate spin system is a mixture of a 
state with the local applied field of magnitude hp = 

The fields are both 



±y/{D + E)/2Eh, for smaU fields, 
parallel and anti-parallel to the distant applied field h, 
but see |J]. (With an equivalent claim for ht which in- 
volves y^{D — E)/2Eht.) Experimentally the existence 
of this superposition is easy to establish. The transverse 
field h couples to the physical 5*2, i.e., this operator in 
the non-singular gauge, and which has a spectrum with 
a period h or in units of magnetic moment gfj-sfi where 
g is the g-factor. Thus if M^ is the measured moment 
in the direction of the transverse field, m = Mz/gfisfi, 
the dimensionless magnetization has a period of unity. 
This would be the case if the state did not consist pre- 
cisely of a superposition of field directions. Because 
of the admixture, the measurable period is reduced to 

Am = nil = \ -fjT^ < 1 which is confirmed by exact 



results but remains to be tested experimentally. 
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